Quantum decoupling transition in a one-dimensional Feshbach-resonant superfluid 
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We study a one-dimensional gas of fermionic atoms interacting via an s-wave molecular Fesh- 
bach resonance. At low energies the system is characterized by two Josephson-coupled Luttinger 
liquids, corresponding to paired atomic and molecular superfluids. We show that, in contrast to 
higher dimensions, the system exhibits a quantum phase transition from a phase in which the two 
superfluids are locked together to one in which, at low energies, quantum fluctuations suppress the 
Feshbach resonance (Josephson) coupling, effectively decoupling the molecular and atomic superflu- 
ids. Experimental signatures of this quantum transition include the appearance of an out-of-phase 
gapless mode (in addition to the standard gapless in-phase mode) in the spectrum of the decoupled 
superfluid phase and a discontinuous change in the molecular momentum distribution function. 



o 

CD 



CO 



I 

C 

O 

o 



oo 

in 
o 
in 
o 



-a 
c 

o 
o 



X 



Recent experimental advances [l| have led to a realiza- 
tion of paired superfluidity in degenerate atomic gases. It 
is driven by the atomic Feshbach resonance (FR) through 
a molecular state whose rest energy (detuning) v can be 
adjusted with a magnetic field. The associated high tun- 
ability of interactions allows one to explore superfluidity 
in these systems ranging from the BCS regime of strongly 
overlapping Cooper pairs (for large positive detuning) to 
the BEC regime of dilute Bose-condensed molecules (for 
negative detuning) |2|. 

In all BEC-BCS crossover studies to date, it has been 
tacitly (correctly 0, 0]) assumed that the superfluid 
phases of the closed-channel Bose-condensed molecules 
and open channel Cooper-pairs are locked together by 
the FR coupling, with the superfluid at low-energies 
characterized by a single gapless Bogoliubov (Oth sound) 
mode. In this Letter we show that in striking contrast, 
a one-dimensional (Id) resonantly interacting atomic 
gas 0, ElEL realized by a sufficiently high aspect ra- 
tio trap [a, 111, can exhibit a more interesting possibil- 
ity. Namely, for a range of parameters quantum fluctua- 
tions, enhanced by the low dimensionality, suppress the 
FR coupling thereby leading to a quantum phase tran- 
sition into a superfluid state where the Cooper-pair and 
molecular superfluids are decoupled. Our main findings 
are summarized by the phase diagram in Fig. ^ One 
striking experimental signature of this decoupling tran- 
sition is the appearance of an out-of-phase gapless mode 
(in addition to the abovementioned standard gapless in- 
phase mode) in the spectrum of the decoupled superfluid 
phase, that should be observable through Bragg spec- 
troscopy ^j. Another signature is a jump across the 
transition in the exponent a characterizing the molecu- 
lar momentum distribution function nf,(k) oc k~ a with 
adecomjied < ^coupled, measurable via time-of-flight im- 
ages |||- 

The decoupling transition can be seen through the 
bosonization representation of the molecular and atomic 
quantum fluids as acoustic charge and spin collective 
modes, with the FR interaction reducing to a Josephson 
coupling between phases of the molecular and Cooper- 



pair superfluids. Now, a sufficiently strong atomic re- 
ulsion (attainable through optical lattice engineering^, 
Hll|). that "localizes" atom number can be tuned to 
enhance the canonically-conjugate superfluid phase fluc- 
tuations, to the point that the FR (Josephson) coupling 
is averaged away at low energies. 




Coupled Spin-gapped (CSG) Decoupled Spin-gapped (DSG) 



FIG. 1: (Color online) Main: Phase diagram as a function of 
the scaling dimension Df T of the FR coupling and the direct 
atomic interaction strength gt s for g{ r = 0.19. For attractive 
interactions, g^ s < 0, fermion spins are gapped and the tran- 
sition is between coupled and decoupled spin-gapped paired 
superfluids. For repulsive atomic interactions, gt s > 0, the 
transition is between the coupled spin-gapped state and a de- 
coupled spin-gapless state. Insets a,b,c: Plots of spectrum 
for fermion spin (s), fermion charge (c), and boson charge (b) 
excitations within these three phases. 

This transition has strong connections to other 
interesting examples of fluctuation-driven decoupling 
of Josephson-coupled XY-models, that fall into the 
roughenning universality class [l^. Most notably, these 
include the sliding phases of DNA-lipid complexes fl^j 
and their quantum higher-dimensional Luttinger liquid 
(LL) generalizations[14|, the latter a realization of a long- 
sought-after LL power-law phenomenology in d > 1. 

We now demonstrate and explore the superfluid de- 
coupling transition through a sketch of our detailed 



2 



calculations 15j. The appropriate model of a two-species 
°~ =1\ I (two distinct values of an internal degree of 
freedom, e.g., hyperfine states) fermionic atomic gas 
interacting through a molecular FR is described by a 
Hamiltonian 2] H = £ fe(T w\„ c h„ +E 9 {e q /2 + y)b\b q + 
Hf r + H m t, where = h 2 k 2 /2m, v is the FR detuning, 
and Hi„± is the direct fermionic and bosonic interactions. 
The FR interaction 

H h = -g iT J d d x(c\(x)c\(x)b(x) + h.c), (1) 

describes the interconversion between a pair of open- 
channel fermionic atoms Cka and a closed-channel di- 
atomic molecule b q . The form of this atom-molecule cou- 
pling guarantees that the Cooper-pair and molecular su- 
perfmids must condense together, and that, in d > 1, 
tuning through an s-wave FR leads to a smooth BEC- 
BCS crossover (rather than a phase transition) between 
two limits of a single paired superfluid state. 

However, in Id long-range order (BEC) is precluded 
even at T = by enhanced quantum fluctuations (lif . 
and superfluidity is characterized by a quasi-long-ranged 
superfluid order, that admits a richer set of possibilities. 
To see this, it is convenient to use a bosonized represen- 
tation [171 Il8| of the atomic and molecular fields, with 



b(x) 



J<t>b(x) 



(2) 



c a (x) ~ e ik 1 x e i(4 >a (x)+9 < ,(x)) _|_ e -ik I x e i(^ a (x)-e <T {x)) ^ ^ 

an approximate form sufficient here. In the above kf is 
the Fermi wavevector, the two terms in Eq. are the 
right and left moving contributions to c CT , and 9 a and 
4> a are the fermion charge-density wave (phonon) and 
superfluid phase fields, respectively. These can be conve- 
niently split into charge (c) and spin (s) components via 
4>c. s — (</>| i 0|)/\/2 and similarly for 9 CS . 

Using Eas. (|2l3|) inside H and focussing on densities 
incommensurate with the optical lattice, we find the 
coherent-state imaginary-time action (with Ti = 1) 

S = J dxdr( ^[Mdxtaf+V^idr^) 2 ] 



+^ rir [v s (d x e s ) 2 + v-\d T e s ) 2 } + 
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9h 



9bs 
27ra 2 



cos 



V89 s 



cos((j)b — V2<pc) cos a/20; 



(4) 



governing the dynamics of three acoustic modes (</>b, <fi c , 
and 9 S ) coupled by two key nonlincaritics gb s and gi T , 
that arise from the backscattering part of fermionic short- 
range interaction gi, s n^(x)ni(x) and from the FR cou- 
pling, respectively. Here, the dimensionless couplings 
9bs = (jbs/nvf and g lr = 2gf T i/p b ~/irvfPf 1 with p f and 
Pb the dctuning-dependent fermion and boson densi- 
ties, Vf — Trpf/2m the Fermi velocity. All other rel- 
evant atomic and molecular interactions have been in- 
corporated into the Luttinger parameters Kb, K c and 



K s , and velocities Vb, v c and v s . For weak interac- 
tions K StC sa 1 ± (?6 S /2 and Kb — > oo. In the above 
a = {t^v 1 / 2 Pf)^ 1 is the ultra-violet (UV) cutoff. 

Based on experience with the sine-Gordon model 0, 
Hif we expect and find that over a large range of Luttinger 
parameters K a a pertubative treatment of the backscat- 
tering and FR nonlinearities breaks down. A renormal- 
ization group (RG) treatment is therefore necessary to 
ascertain the low-energy behavior of the system. This 
amounts to successively integrating out the high-energy 
degrees of freedom in an infinitesimal wavevector shell 
a~ x e~ < q < a -1 around the UV cutoff. Resulting ef- 
fective couplings appearing in S then determine the low- 
energy thermodynamics of the system. 

Before discussing the full behavior it is convenient to 
consider a limiting regime of the above model. The sim- 
plest case is that of large attractive atomic interactions, 
gbs -C — 1 (that, as we will see below, is induced by a 
finite FR scattering even if the nonresonant atomic in- 
teraction gbs is moderately repulsive). In this limit, the 
backscattering nonlinearity "freezes" 9 S at 0, correspond- 
ing to a spin-gap state with gap w 2y/\gb s \v s K s /a. 
Using 9 S = inside S reduces the problem to two acoustic 
modes <t>b, c coupled by a single FR nonlinearity. Its effect 
can be assessed using an RG analysis that here amounts 
to a determination of the scaling dimension D of the 
FR operator 0(x) = cos((f>b(x) — \f2(j) c (x)) cos V29 s (x) 
around the Gaussian fixed point, gb s = 0,<7f r = 0. 
This in turn is determined by the long-scale behavior 
of (O(x)O^(0)) ~ \x\- 2D . A simple calculation in the 
spin-gap phase gives D = Df T = (4Kb)~ 1 + (2if c ) -1 . A 
standard RG analysis then shows that for Df r < 2 
(Df T > 2) the effective FR coupling grows (diminishes) 
at low energies under the RG coarse-graining procedure. 
For D[ v < 2 and weak <?f r <C 1, on scales longer than 



ps aVf g fr tr , the FR coupling dominates over 
the kinetic energy and the growth of gt r saturates at 

2 

g^. b . In this strong-coupling regime, the FR interac- 
tion locks the closed-channel molecular superfluid phase 
4>b to the open-channel Cooper-pair superfluid phase 
V2(f> c , a characteristic of the "coupled spin-gap" (CSG) 
state. Approximating the FR coupling by a harmonic 

2 

spring {gf r D{r /2ira 2 )((f>b — V2(j)c) 2 , with a "stiffness" soft- 
ened by quantum fluctuations up to the scale £f r , and 
diagonalizing the quadratic form in <pb and <f) c , gives two 
dispersions of the CSG state: 

1 



A 2 

^fr 



T^J[{v 2 +v 2 )k 2 + 



A 2 I 2 



Av%v 2 k 4 



4c 2 PA f 2 r 



(5) 



with Afr = a 1 \gir\ 2 ~ D ft- y VbK b 1 + 2v c K c 1 , and c = 

\J (v c K c + 2v b K b )/(vc 1 K c + 2v~ b 1 K b ). The dispersion 
w+(fc) (~ ck) characterizes the gapless Bogoliubov mode 
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corresponding to in-phase oscillations of the closed- 
channel molecular (ip b ) and the open-channel Cooper-pair 
(4> c ) superfluid phases. The dispersion w_ (« A fr ) is for 
the gapped mode (in the long wavelength limit given by 
= <t>b — V2<t>c) in which (f> b and (f> c fluctuate out phase. 
Hence, the CSG state is characterized by two gapped 
modes 6 S , <p- and a single in-phase gapless Bogoliubov 
mode (see Fig. \%)p) . In this coupled state the molecular 
phase 4>b is characterized by an effective Luttinger param- 
eter K b = yj {v b K b + v c K c /2)(v- 1 K b + Vc 1 K c /2) > K b 
that can be read off from the action, Eq. (0J after 
simply imposing the low-energy FR coupling constraint 
<i>c = 4>b/V2 inside S. Similarly, atomic charge correla- 
tions are controlled by K c = 2K b . 

Now consider Df r > 2 but still deep within the spin- 
gapped state, gbs <C —1. Here, quantum fluctuations of 
4>b, c become sufficiently strong so as to average away long- 
scale effects of the FR coupling, reducing it relative to the 
kinetic energy of fluctuations on these long scales. We re- 
fer to this distinct thermodynamic state (special to Id) 
as the 11 decoupled spin-gap" (DSG) superfluid. It is char- 
acterized by effectively independent fluctuations of the 
low energy molecular and Cooper-pair superfluid phases 
and therefore exhibits two gapless superfluid modes dis- 
played in Fig. ^ (with velocities v b and v c ), observable 
via Bragg spectroscopy [Tof . 

Hence through the decoupling transition charge Lut- 
tinger parameters "jump" from K biC down to K bjC . This 
is reflected in the corresponding momentum distribution 
functions n b .c(k), that for molecules [using Eq. is 
given by n b (k) oc / dxe* kx (M '(x)b(0)) ~ k~ a (for k -> 0), 
with a C oupicd = 1 - {2K b y 1 > a d ccoupiod = 1 - (2K b y 1 . 
This abrupt enhancement in the low k peak of Ub{k) re- 
flects the suppression of molecular phase fluctuations in 
the coupled phase (due to locking to 6 C ) and should be 
measurable via time-of- flight images |8| • 

Hence, as advertised in Fig.l, we predict a quantum 
phase transition between the coupled, CSG and decou- 
pled, DSG spin-gapped superfluids. To determine the 
phase diagram outside of the deep spin-gap regime, i.e., 
for values of the backscattering amplitude gb s other than 
gbs *C — 1, requires a detailed RG analysis. As outlined 
above the RG computation involves progressively inte- 
grating out (perturbatively in g^ and g bs ) degrees of free- 
dom at the UV cutoff scale a -1 . For simplicity we spe- 
cialize here to the case of equal velocities (v b ~ v c ~ v s ), 
leaving the more technically challenging general case to 
a future study [lij . The result of this coarse-graining 
procedure is summarized by the RG flow equations 



dgbs 
di 

dgh 
dl 



2g hs {\ - K s 



\9%: 



(2 - Ar - Mg^ 



(6) 
(7) 



dl 
dDfr 
~df 



1 2 1 2 
"~j9bs ~ 4#fr) 

9 2 

■gfc. 



(8) 
(9) 



One important subtlety |l8j is that abelian bosonization 
does not explicitly capture the underlying SU(2) spin- 
rotation symmetry embodied in our system's Hamilto- 
nian. It can, however, be restored (as we have done 
above) by a specific choice of the UV cutoff a and by 
imposing the relation K s = 1 + g bs /2 on the initial con- 
ditions for the flows. This constraint stems from the 
fact that the value of g bs and the correction to K s arise 
from the same fermion interaction. It is simple to ver- 
ify that this SU(2) invariance relation between K s and 
gbs is preserved by the RG flows Eq. JBJ and ©, i.e., 
that for K s = 1 + g bs /2 initially, d(l + gbs/2 - K s )/dt = 
g bs (l + g bs /2-K s ) = 0. 

Incorporating this relation significantly simplifies the 
RG flow equations, giving: 



dgbs 
di 

dgh 
di 



-9bs 



5f r /2, 



(3/2-D ft -3 5fcs /4)<7fr, 



(10) 
(11) 



along with Eq. © which is unchanged. The phase dia- 
gram Fig.^is determined by the asymptotic (large I, cor- 
responding to low energies and long scales) flows of gh-(i), 
D( T (l) and gbs(i)- A priori, one might have expected four 
distinct phases corresponding to four different combina- 
tions of relevent and irrelevant regimes of two couplings 
g bs and gf T . However, as is clear from the structure of 
the flow equations [particularly Eq. (fTUf) ] . it is impossible 
to have a phase-coupled but spin-gapless state charac- 
terized by an asymptotically nonzero gh{i — > oo) and a 
vanishing g bs . That is, back-scattering g bs (that leads to 
Cooper-pair singlet formation) is always relevant when 
gtr is. Physically this can be understood by observing 
directly from S that an arbitrary strength FR coupling, 
at T — always induces pairing in an itinerant (i.e., ig- 
noring Mott-insulating effects of a commensurate lattice) 
fcrmionic atom system. 

One simple limit of the flow equations is D{ r 3> 1, in 
which case gf r clearly flows to and the phase boundary 
separating the DSG (for g bs < 0) and decoupled spin- 
gapless DNG (for gbs > 0) states asymptotes to g bs = 
0, respectively corresponding to attractive and repulsive 
nonresonant atomic interactions. 

Another observation is that, for large repulsive back- 
ground interactions g bs 3> 1, the phase boundary sep- 
arating the CSG and DNG states asymptotes to Df T = 
3/2, as can be seen from the flow equations. These show 
that a strong background repulsion strongly suppresses 
the FR coupling g^ through Eq. (fTT|) , so that it has little 
effect on the flow of g bs , which then itself can flow to 
along the gf T ~ sine-Gordon separatrix. This then re- 
duces the eigenvalue for gf r in Eq. (|ll|l to 3/2— Df r leading 
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to the phase boundary Df r = 3/2. This is consistent with 
the observation that the FR coupling scaling dimension 
D = D{ T + \K S is 2 at the transition and K s = 1 in the 
DNG state, constrained by the spin-rotational invariance. 
For an intermediate and small positive values of gbs , we 
have numerically integrated the RG equations to deter- 
mine the interpolation of the vertical part (£>f r = 3/2) to 
the horizontal part (gbs = 0) of the phase boundary, as 
illustrated in Fig. ^ 

In contrast to repulsive background interactions that 
can be irrelevant if the FR coupling is, attractive interac- 
tion (gi, s < 0) always grows under coarse-graining. Physi- 
cally, this is a reflection that Cooper-pairing always takes 
place at T = 0, even in the absence of the FR coupling 
that can only enhance spin-singlet formation. Hence, at 
sufficiently long scales, the RG flows leave their perturba- 
tive regime of validity confined to \gb s \ < 1 and |<?f r | < 1. 
Outside of this range requires a nonperturbative analysis. 
As discussed earlier, the deep spin-gap state gbs <C — 1 
can nevertheless be analyzed by simply setting 8 S = 
in the action and recomputing the flows perturbatively 
in gf r . As illustrated in Fig. ^ this leads to a phase 
boundary between CSG and DSG that is asymptotically 
vertical and given by D[ r = 2. 

To calculate this phase boundary for more moderate 
(but still negative) values of gb s requires a "matching" 
calculation. To see this we note that for small nega- 
tive gbs, perturbative RG flows (describing spin fluctu- 
ations on shorter length scales inside the spin-gapped 
state), remain valid up to a crossover spin-gap length 

scale, £bs ~ avl^ 2 \gbs\~ 1 ^ 2 ~ 2K ^ (coinciding with the 
width of the soliton in the sine-Gordon model). Hence 
we integrate the perturbative RG flows out to length 
£bs = avs e** , defined by gb s (i*) ~ 1, beyond which 
spin fluctuations freeze out, 9 S ~ 0. On longer scales we 
match onto the strongly coupled spin-gapped state, char- 
acterized by setting 8 S = inside the action, Eq. Q but 
with the fiuctuation-renormalized parameters gh(£*) and 
D b (t) determined by the RG flow. Setting D b \t*) = 2 
then determines the fiuctuation-renormalized CSG-DSG 
phase boundary illustrated in Fig. ^ The shift of the 
phase boundary toward a larger critical Df v can be seen 
from the decrease of D{ T (£) under coarse-graining, which 
in turn corresponds to enhancement of superfluidity by 
the FR coupling. 

A nontrivial challenge is the experimental realization of 
a LL with parameters tunable across the decoupling tran- 
sition. The necessary large value of Df r > 3/2, requires 
small Luttinger parameters, K a , realized by strongly in- 
teracting systems, for which a relation to microscopic 
parameters is difficult to establish. One exception is the 
extended quarter-filled Id Hubbard model, where in the 
limit of large on-site repulsion, K c = [2 + A/it sin -1 v]^ 1 , 



with v = V/(2\t\) a ratio (limited to < 1 by phase sep- 
aration) of the nearest neighbor hopping t and repul- 
sion V energies, respectively jj^- Hence, K c as low 
as a 1/4 can be reached, assuring a transition (that for 
strong repulsion takes place at Z?f r = 3/2) into the de- 
coupled superfluid state at v — 1/V2 for weakly inter- 
acting molecules. In the more favorable molecular Tonks 
regime with Kb — > 1, we predict the quantum transition 
at v « 0.38. Thus we suggest that a Feshbach-resonantly 
interacing atomic gas confined in a highly anisotropic 
(Id) trap and subjected to a periodic optical potential @ 
is a promising candidate for an experimental realization 
of the phase diagram and the decoupling transition dis- 
cussed here. 
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